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I. INTRODUCTION 
According to Haldar and Mahadevan 


uncertainty. In general, repeated measurements of physical 
phenomena yield multiple results. The occurrence of 
various outcomes with no pattern is described by terms 


such as uncertainty and randomness. 


According to Real et al. [2], the main quantities 
involved in engineering design are, in reality, random 
variables, which have a probability distribution. Therefore, 
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[1], 


observable phenomena contain a certain amount of 


Abstract — Safety is an essential requirement of a structural system. 
Reliability is an additional tool of growing importance in engineering, as it 
allows us to quantify uncertainties in the design. Thus, reliability assists us 
in making more suitable decisions regarding the safety of a structure. The 
present work compares and analyzes structural reliability methods applied 
to various examples of limit state functions. These methods are essential 
tools for this analysis because they identify and quantify uncertainties in 
random variables, allowing the evaluation of the probability of failure of 
the structure. Structural reliability methods were programmed and 
simulated in the Python language. The performance of these methods was 
analyzed through examples of linear, nonlinear, implicit, and explicit limit 
state functions. The results indicate that the simulation of Monte de Carlo 
brute force (MCBF) and importance sampling (MCAI) proved to be quite 
efficient for the examples studied in this work, with values equal to or very 
close to the reference values from the literature. The First Order and 
Second Moment Method (FOSM) presented limitations in some examples 
when the basic random variables do not have a normal distribution and the 
limit state function is nonlinear. The first-order reliability method (FORM) 
employs a failure surface linearization, which does not work well for highly 
nonlinear problems. The second-order reliability method (SORM) has 
improved the FORM results by including additional information about the 
curvature of the limit state function. 


the structure's response to loading will also be a random 


ost variable. 


Thus, ensuring adequate safety levels in structures, and 
providing good performance, is one of the most significant 
engineering challenges. Models based on limit states are 
used in the design standards, where a series of safety 
criteria is specified in the calculation process of structures. 
The effect of uncertainties on the safety of structures can 
be diminished by analyzing the probability of failure of a 
limit state. (ALBUQUERQUE AND REAL [3]). 
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For Verma et al. [4], reliability is defined as the 
probability of trouble-free performance under established 
conditions. Still, it can be considered the probability that 
an element can perform its intended function for a 
specified time under pre-established design conditions. 


After Beck [5], due to the uncertainties of the 
parameters of engineering projects, it is necessary to carry 
out reliability analyzes because of possible failures. The 
answer, in this case, is the determination of the failure 
probability or reliability index associated with a limit state 
function that depends on the random variables of the 
problem. That way, structural reliability methods are 
essential tools for this analysis, as they identify and 
quantify uncertainties in random variables. 


According to Santos [6], several methods have been 
developed to estimate the failure probability and reliability 
index values. Among these methods, FOSM (First Order 
Second Moment), FORM (First Order Reliability Method), 
and SORM (Second Order Reliability Method) stand out. 
These structural reliability methods, also known as 
transformation methods, approximate the limit state 
function at the design point by a hyperplane, in the case of 
FOSM and FORM. In the case of SORM, a quadratic 
hypersurface is used. The approximation of these methods 
can lead to errors, especially when the limit state function 
is strongly nonlinear. In this way, Monte Carlo simulation 
can estimate the probability of failure and the correct 
reliability index, as this method is considered the most 
exact. 


The Brazilian Association of Technical Standards 
(ABNT) is the main standardization body in the country. 
The Brazilian Norms (NBR) are technical documents used 
to standardize and guarantee a quality level. There are still 
no technical specifications based on structural reliability 
principles in Brazil, although this is a reality in several 
other countries. Therefore, this work is justified by the 
omission of the NBR concerning structural reliability. This 
paper aims to compare structural reliability methods 
applied to examples of limit state functions and, thus, 
analyze the efficiency of the method application. 


Il. MATERIAL AND METHOD 


The methodology used in the elaboration of this work 
consists of the development of codes for structural 
reliability analysis through the methods Monte Carlo Brute 
Force (MCBF), Monte Carlo with Importance Sampling 
(MCIS), FOSM, FORM with the HLRF algorithms 
(Hasofer, Lind, Rackwitz and Fiessler) and iHLRF 
(improved Hasofer, Lind, Rackwitz and Fiessler) and 
SORM. 
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After programming, a comparative study between the 
methods was presented, where eighteen examples of 
linear, nonlinear, and implicit limit state functions 
obtained by some researchers are considered. Then, the 
performance and shortcomings of each method are 
analyzed. 


2.1 — Limit State 


Beck [5] states that each distinct form of arriving at an 
undesirable state is generically called a failure mode. 
Where each failure mode gives rise to a limit state 
equation g(X), which is a function of the vector of random 


variables X = [NX Xz. ox.» M_}, in such a way that the limit 


state is determined by equation 1: 


gx) = gX —.X,)= 0 (1) 


a limit state equation g(X) is written in such a way as to 


divide the problem domain into failure My and safety 


domains ft,, as shown in equations 2 and 3: 


a = {xlg(x) < 0) 2) 
n= ilga) > 0) 6) 


The failure domain My is constituted of all points in the 


sample space X that lead to the structure's failure. The 
secure domain fl, is the complementary set to the failure 
domain. 

2.2 — Limit state function 

According to Melchers and Beck [8], just as there are 
uncertainties in engineering projects, a structural system 
has the risk of not reaching the designed performance. 
Consequently, the performance of this system cannot be 
guaranteed to be completely secure. Thus, for Santos [6], it 
is convenient to define the limit state function; considering 
a structure with resistance R and load S, being possible to 
evaluate the failure probability. A limit state function in a 
probabilistic analysis can be determined by equation 4: 


g(x) = R(x) — S(x) (4) 


Where x is the vector of random variables of the problem; 
R(x) is the random variable representing the strength of 


the structure; Sæ} is the random variable representing the 
loading. 
The probability of failure (Fp) is defined when the 


structure's resistance is less than the loading, and the safety 
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probability (F4) is when the resistance is greater than or 


equal to the loading, being given by equations 5 and 6: 


P = P[g{x) <0] (5) 
R =1-}. (6) 


In compliance with the JCSS [8], the probability of 
failure must be calculated based on the type of joint 
distribution of the basic variables and on the standardized 
distributive formalism of dealing with model uncertainty 
and statistical uncertainty. 


Figure 1 shows a surface with the failure region, the 
safety region, and the limit state function for a system with 


two variables. 
g(X) =0 
ri Limit state function 








g(X) <0 


g(X) >0 Failure region 


Safe region 


Fig. 1: Limit state function for a system with two 
variables (ALBUQUERQUE [7]) 


For Beck [5], the reliability index calculated in 
problems involving statistically independent normally 
distributed random variables and a linear limit state 
function corresponds to the smallest distance between the 
origin and the limit state function and is given through 
equation 7. 

as ee (7) 
Voto 

According to Melchers and Beck [8], the failure 
probability can be obtained directly from the reliability 
index for linear limit state functions. It can still be defined 
if the function is nonlinear, but only for an approximate 
(hyper)tangent plane. In each case, # represents the 
smallest distance between the origin in standardized 
normal space and the tangent plane, as shown in Figure 2. 
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Fig. 2: Representation of the reliability index 
(ALBUQUERQUE [7]) 


In Figure 2, the minimum distance point on the 
hyperplane is denominated as the design point. This point 
has the highest probability density within the space 
covered by the failure region, according to Melchers and 
Beck [8]. 


2.3 — Reliability Methods 


Monte Carlo simulation is a numerical method that is 
based on random sampling. Second Beck [5], in structural 
engineering, Monte Carlo simulation adapts to linear and 
nonlinear, static or dynamic problems. In structural 
reliability, the simulation solves, with the same ease, 
problems involving a single limit state equation, time- 
dependent or not. 









Define random variables; 


Generate n samples for the sampling 
function; 


Generate N samples for the random 
variables according to the sampling 
function; 


Solve the function g(x) for each set of 
random variables; 


Extract from the samples the number 
of times the function g(x) = 0; 







Determine the reliability index and 
failure probability. 





Fig. 3: Flowchart of the brute force Monte Carlo 
method algorithm (ALBUQUERQUE [7]) 
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The brute force Monte Carlo simulation method 
(MCBF) is a numerical technique based on repetition of 
the simulation process, using a specific set of values in 
each simulation. 


According to Haldar and Mahadevan [1], the MCBF 
simulation has six essential elements, and the algorithm of 
the method is presented in Figure 3. 


One of the techniques to accelerate the convergence 
of Monte Carlo simulation is utilizing importance 
sampling. The technique determines the probability of 
failure in the limit state region. Beck [5] states that the 
importance sampling technique uses a sampling function 
to avoid generating samples far from the area of interest, 
far from the failure domain. The points are generated from 
a sampling function hy(x). Multiplying and dividing 
Equation 5 by hy(x), we get equation 8: 

1%” yf) (8) 
Bx R=—) Lx) = 
f f N = o (x; hy læ) 
The point x; is associated with a weight wẹ == 1, 
represented by the ratio of the functions f(x) and hy(x). 


The variance of the failure probability is estimated by 
equation 9: 


E 1 N a 
vertõi= qr) Day olsen A) 0) 


According to Haldar and Mahadevan [1], the algorithm 
of the MCIS method is presented in Figure 4. 


Define random variables; 


Use the design points obtained 
through the FORM iHLRF method; 
















Generate n samples for the sampling 
function; 


check the occurrence of failure for 
each sample, through the indicator 
nction; 










Estimate the average failure 
probability through Equation 8; 


Estimate the variance of the failure 
probability through Equation 9; 


Determine the reliability index and 


failure probability 





Fig. 4: Flowchart of the Monte Carlo method 
algorithm with sampling by importance (ALBUQUERQUE 
[7]) 
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According to Haldar and Mahadevan [1], FOSM and 
FORM can be used to evaluate a function of normal and 
statistically independent variables. FORM can also 
represent a nonlinear limit state function by the first-order 
approximation with equivalent normal variables. And 
SORM determines the probability of failure by 
approximating the nonlinear limit state function by a 
second-order representation. FORM substitutes the limit 
state function at the design point by a hyperplane, while 
SORM uses a quadratic hypersurface. In Figure 5, the first 
and second-order approximations are shown. 





Fig. 5: First and second-order approximation - FORM 
e SORM (ALBUQUERQUE [7]) 


According to Beck [5], in the FOSM and FORM 
methods, the limit state equation is approximated by a 
linear function. The information for constructing the joint 


probability density function f(x) is limited to moments 


up to second-order (mean and covariance). A 
representation of the problem variables by their moments 
corresponds to assuming them as having a normal 
distribution. 


The FOSM method solution algorithm is presented 
below, following (BECK [5]). 


Algorithm 1. Reliability method FOSM 


1. choice of initial failure point x™ to k=0 


(usually the mean point); 
2. evaluation of Jacobian matrices (yx and 


Ixy) and the vector of the means V; 
point transformation x* of E — ¥; 


evaluation of g (x*); 


So ae SM 


gradient calculation: 
a. calculating the partial derivatives of g(x) 
in design space E 
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b. transform the gradient to ¥ 
c. calculating director cosines a (y*) 
6. calculation of the reliability index (2) 


7. calculation of the new point of failure 


y** through 


yi? = phi — amos, (10) 


8. transformation of O ago in E; 


9. verification of the convergence criterion. 
If: 


Vely***) ye 
1————_—_ DT 
live] [ly] 


ul 
<e ekG™)i<s O 


being e = 10~*; the algorithm is interrupted (y* = 
yet ); otherwise, return to item 4 with k=k+1 


until reaching convergence; 
10. at the end, evaluation of the reliability 
index at the design point: B = || y*|). 

The FORM algorithm is equal to the FOSM for a 
problem with a normal probability distribution of random 
variables. When distributions are mixed, it is necessary to 
calculate the mean and standard deviation of an equivalent 
normal distribution (ALBUQUERQUE [7]). 






g(y)=c>0 


Fig. 6: Representation of the FORM method with 
the HLRF algorithm - Linear limit state function 
(ALBUQUERQUE [7]) 
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According to Melchers and Beck [8], the search for the 
design point y* is one of the fundamental steps to 
obtaining the failure probability by the FORM method. 
The transformation of any variables, correlated or not, into 
statistically independent equivalent normal variables is 
performed through the Nataf model transformation. 

For Beck [5], the solution of reliability problems via 
FORM for nonlinear limit state functions involves 
resolving an optimization problem to search for the design 
point. The most used algorithm is Hasofer, Lind, Rackwitz 
and Fiessler (HLRF). For nonlinear limit states, the 
minimum distance calculation becomes an optimization 
problem. The solution to this problem is presented in 
equation 12: 


a= jlyl] = V¥ y? (12) 


Figures 6 and 7 illustrate an iteration of the FORM 
method with the HLRF algorithm, considering a linear and 
nonlinear limit state function. It is observed that while the 
design point is not reached, the condition gly) = 0 is not 
obtained. That is, the point y* does not reach the failure 
surface, but a surface to which gly) is constant (c). 
Therefore, the optimization algorithm must start from the 
point y*, which may not be in the limit state, and converge 


to the minimum distance point y* in the limit state. 





Fig. 7: Representation of the FORM method with the 
HLRF algorithm - Nonlinear limit state function 
(ALBUQUERQUE [7]) 


k+l can be obtained 


The recurrence formula for y 
considering the linear approximation of gy), and is 


represented through equation 13: 
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[PgG*)'y* — gy FP gG*) (13) 
rg GI " 
According to Sudret and Der Kiureghian [10], the 
HLRF algorithm can be improved (IHLRF) through a 


linear search for the fit. The original HLRF algorithm 
gives the search direction of equation 14: 


_ WgG*¥ y* — g0") 
lira 


The solution algorithm of the FORM iHLRF method is 
presented following (BECK, [5]). 


Algoritmo 2. Reliability method FORM iHLRF 
Data: y¥ e R”,A > On > La e 0.10 < m,<1 


pert = 


(14) 


dë = yet ye Vg (y*) ae y* 


1. Calculate the search direction of dE 


given in equation 14; 
2. Determine the penalty parameter c of the 
merit function 


fley") > a 








c=nmax{ PI silva" 
1 les Gr] Teto] 
If not 
con W. 
7 te 


3. Determine cf using inequality 


While p(y* + ad*) — p(y®) > m,a¥p(y*)' at 


Do a= 


NIA 


Č = a; 
4. Update the point y: 


The SORM approximates the limit state function at the 
design point by a quadratic hypersurface. Beck [5] states 
that this method results in a good approximation of the 
failure probability when the limit state function, in 
standard normal space, is linear or weakly nonlinear in the 
vicinity of the design point. 


For Cordeiro [11], the SORM method estimates the 
failure probability through a second-order approximation 
of the limit state function at the design point. The second- 
order derivative of this function approximates the 
curvature of the limit state function with respect to the 
equivalent normal space variables. The SORM method 
algorithm is calculated by obtaining the reliability index 
and failure probability of the FORM method. The 
algorithm of the FORM and SORM methods is shown in 
Figure 8. 
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According to Haldar and Mahadevan [1], the curvature 
of any function is related to the second-order derivatives in 
relation to random variables. Therefore, the second-order 
reliability method (SORM) improves the FORM result by 
including additional information about the curvature of the 
limit state of the function. 


HI. RESULTS AND DISCUSSION 


Eighteen problems are presented for comparative 
analysis, applying Monte Carlo simulation and 
transformation methods in structural reliability problems. 
The study consists of a comparative analysis of the 
convergence of the reliability index of linear, nonlinear, 
and implicit limit state functions found in the literature. 


The transformation methods are interrupted when they 
encounter a point whose stationarity measure value is 
sufficiently small. In this way, the conditions are 
established according to equation 11, using E= 10€. 
The method execution is interrupted if these conditions are 
not met until the code reaches 100 iterations. 

In implementing Monte Carlo simulation, MCBF 
works well for a simulation number up to 9x 107 . The 
MCIS uses a predetermined point as sampling, using the 
design point obtained through the FORM iHLRF method. 

The probability distributions, the moments (mean and 
standard deviation) of the random variables, and the limit 
state function gix) are presented for each problem 
analyzed 

Problem 1 - Beck [5]: Limit state function 
glx) = x,x,— xy, the variables follow normal probability 
distribution, with mean values uy = 40, uy, = 50 and 
Hx, = 1000 and standard deviations Gy, =5,0y, =2.5 
and Sy, = 200. 

Problem 2 - Shayanfar [12]: Limit state function 
glx) =x? + xix, + x3 — 18, the variables follow normal 
probability distribution, with parameters uy = 10, 
By, = 9.9 and oy, =oy, = 5. 

Problem 3 - Keshtegar e Chakraborty [13]: Limit state 
function glx) =x? + x — 50, the random variables, x, 
has lognormal probability distribution with parameters 
Ux, = 5 and Sy, = L, xz follows Gumbel distribution with 


parameters Hy, = 10 and o, = 10. 
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(start ) Define the limit state function g(x): 


Choose a point of failure x” initial; 


+ 
Assume an initial value of the reliability index (8); 


























Point transformation x” de X — Y; 





x . e r ð, 
Calculate Partial Derivatives E 
i 


Calculate director cosines («,,) at the design point; 








Determine the new failure point with: y’ = Hri = ay Bok, 


Verification of the convergence criterion; 








Calculate the P Frome e Ê Form- 


Orthogonal matrix calculation; 
z 
Determination of second order partial derivatives E 
2 
y, 






Calculation of curvatures in the region of the most likely point of failure; 





Calculate the Pf ogm € Bsorm- (ne ena) 


Fig. 8: Algorithm of the FORM and SORM methods (ALBUQUERQUE [7]) 


Problem 4 - Li e Zhang [14]: Limit state function 
g(A,P) = 30 — u,(A, P), the random variables, x, has 


lognormal probability distribution with parameters 
Hy, = 0.0025 and o, = 0.00025, x, follows Gumbel 


distribution with parameters Hy, = 200 and gy, = 50. The 





authors present an analysis of the plane truss with the 
dimensions and loads shown in Figure 9. The bars' 


material has a linear elastic behavior, with a modulus of Fig. 9: Structural scheme of the flat truss (Li e Zhang 


elasticity E = 200 GPa. The cross-sections of the truss H41) 

members have an area of 2Z5cm?. The structure is 

subjected to loads P = 200 kN, applied vertically at nodes Problem 5 [Ang e Tang (2007)]: Limit state function 
2 and 4. g(f.S) =R —G—Q—W, the variables have a normal 


probability distribution, with parameters 
Ha = 975, we = 200, ug = 300, uy = 150, op = 
146.25, og = 14, og = 36 and oy = 30. 
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Problem 6 - Ang e Tang [15]: In this example, the 
same data from Problem 5 are used, with changes in the 
distributions of the variables, R and G have lognormal and 
normal probability distributions, respectively, Q and W 
follow the Gumbel probability distribution. 


Problem 7 - GHALEHNOVI et al. [16]: Limit state 
function g(x) =exp(0.2x, + 1.4)— x, the variables 
follow normal probability distribution, with parameters 
Hy, = By, = Dando, = Gy, = 1. This function is weakly 


nonlinear and its behavior is illustrated in Figure 10. 
7 


01 



















Fig. 10: Failure threshold and failure and safety 
domains, for problem 7 (ALBUQUERQUE [7]) 


Problem 8 - SANTOS [6]: Limit state function 
g(p.q) = 3— q + (4p)*, the variables follow normal 


probability distribution, with parameters gy = uy =0 and 
Sp = @, =1. This function is strongly nonlinear and its 


behavior is illustrated in Figure 11. 





| | 
— 9(p,q) = 0 
33 


3.0 + 














Fig. 11: Failure threshold and the failure and safety 
domains, for problem 8 (ALBUQUERQUE [7]) 


Problem 9 - Beck [5]: Limit state function 
g(x) =x, +2x,+3x,+x,—Sx5x,, the 


follow lognormal probability distribution, with mean 


variables 
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values ty, = Hy, = Hy, = By, = 120, py. = 50 and 


Ex. = 40 and standard deviations 
to 


Tx, = Oy, = Oy, = Gy, = Gy, = 12 and oy. = 15. 


Problem 10 — GROOTEMAN [17]: Limit state 
function g(x) = 25-— 0.2357(x,—x,) +0.0046(x, + x, — 20)º, 
the variables have a normal probability distribution, with 
parameters Hy, = Hy, = 10 and o,, = Gy, = 3. 

Problem 11 - SANTOS et al. [18]: Limit state function 
g(x) = x,x, — 78,12x,, the random variables, x, and xz 
have lognormal probability distribution with parameters 
Hy, = 2x 10”, uy, = 10º, o,, = 0.5 x 10º and 
Ty, =0.2x 10-74, x, follows Gumbel distribution with 
parameters gy. = 4 and Sy, =1. 

Problem 12 - SANTOS et al. [18]: Limit state function 
glx) =xi+xi— 18, the variables follow normal 
probability distribution, with parameters gy = gy, = 10 
and Gy, = dy, = 5. 

Problem 13 SANTOS et al. [18]: Limit state function 
g(x) =x] +x3—18, the variables have a normal 


probability distribution, with mean values py, = LO and 
Hy, = 9,9, and standard deviations Gy, = Gy, = 5. 

Problem 14 - Choi et al. [19]: Limit state function 
g(x) =x}+2x3—20, the variables have a normal 
probability distribution, with parameters gy = Hy, = 10 
and gy, = Gy, = 5. 

Problem 15 - GHALEHNOVI et al. [16]: Limit state 
function glx) = =x 1º + 2x3 + x, +3, the variables have 
a normal probability distribution, with parameters 
Hy, = Hx, = Hx, = O and ay, = Gy, = Gy, = 1. 


Problem 16 - [KESHTEGAR e MENG, 2017]: Limit 
state function 


gx) = 0.03 — (h) (g) +0) 


the random variables, q.LA, and À. have normal 
probability distribution with parameters u; = 20000, 
S = 1400, m=12, 9 =0.12, ua =982x 10, 
Sa, = 3.9852 x 10-5, Ha = 0.04, qq = 00048, E, 
follows lognormal distribution with parameters, 
Be, = LO x 10“, oz, = 60x 10º and E. has Gumbel 


distribution with parameters, uz =2.0x 10!º and 
Sg = 12x 10°. In this example, a roof truss under 


uniform load is analyzed, shown in Figure 12. The upper 
chords and compression bars of the truss are made of 
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reinforced concrete, and the bottom chords and tension 
bars are made of steel. 





| 0.2781 | 0.2221 | 0.251 | 0.251 | 


Fig. 12: Vista esquemática da treliça de telhado 
(KESHTEGAR e MENG [20]) 


Problem 17 - [SHAYANFAR et al., 2017]: Limit state 
function g(A,. A, 4,.B,P,E)J=01-B. u(A, A, dy, P,E), 
the random variables, Apn eB have normal 
probability distribution with parameters ma, = 107, 
Ga, =5.0x10*, pa =1L5x107, o, =7.5x 107, 
Ha, = 60x107, o, =3.0x10~, pg = 1.0, og =0.1, 
P follows lognormal distribution with parameters, 
up =25.0x 10, g,=25x10º and E has Gumbel 
distribution with parameters, up=69x10* and 
ge = 3.45 x 107. In this problem, a 10-bar truss structure 
is presented in Figure 13 with three different cross- 
sectional areas A, (vertical), A, (horizontal) e As 
(diagonal). The structure is subjected to two external 
horizontal loads P. The failure of the structure is 
considered when the horizontal displacement in the upper 
right comer exceeds the allowable displacement 
d,= 0.1m. is the modulus of elasticity. The random 
variable B was introduced for the model uncertainties and 
L = 9.0 m is a deterministic variable. it is assumed that Ay, 
A, , Ay, B, P, E are statistically independent random 


variables. 
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Fig. 13: 10-bar truss structure (SHAYANFAR et al. 
[12]) 


Problem 18 - ROUDAK ef al. [21]: Limit state 

‘ x 
function glx) = Cso) — 360 ( cad the random 
variables, x4, Xz e xy follow normal probability 
distribution with parameters uy = 10.0, q, = 0.4, 
Hy, = 20x 10’, q, =0.5 x10’, 


Sx, = 0.01. In this problem, a continuous beam with three 


Ex, = 0.4 e 


spans is studied, considering the parameters uniform load 
intensity, modulus of elasticity, and section height. Figure 
14 illustrates the analyzed beam. 





7m 7m 7m 


Fig. 14: Three-span continuous beam (ROUDAK et al. 
[21]) 


Table 1 shows the results of the reliability indices of 
each method studied for each problem. 
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Table 1 - Results of the reliability indices of the analyzed problems 






























































Methods 
Problems MCBF MCIS FOSM FORM FORM SORM 
HLRF iHLRF 
1 3.04 3.04 3.04 3.04 3.04 3.04 
2 2.52 2.52 1.54 1.54 2.30 2.62 
3 3.56 3.56 1.51 1.84 3.26 3.53 
4 3.20 3.20 4.12 3.19 3.19 3.19 
5 2.10 2.10 2.10 2.10 2.10 2.10 
6 2.33 2.35 2.10 2.45 2.45 2.45 
7 3.38 3.38 3.35 3.35 3.35 3.38 
8 3.57 3.57 3.00 3.00 3.00 3.00 
9 2.94 2.94 3.65 3.04 3.04 3.04 
10 2.75 2.75 2.50 2.50 2.50 2.50 
11 4.43 4.43 3.29 4.43 4.43 4.43 
12 2.55 2.55 2.24 2.24 2.24 2.63 
13 2.55 2.55 1.16 1.16 2.22 2.63 
14 2.92 2.92 0.93 0.93 2.36 2.84 
15 3.43 3.43 2.52 2.52 3.00 3.37 
16 2.35 2.35 0.18 0.59 2.42 2.42 
17 4.30 4.30 5.49 3.99 4.30 4.30 
18 2.52 2.52 not not 2.52 2.52 
converged | converged 























According to Table 1, problems 1 and 5 tend to 
converge to the same point, represented by the reliability 
index. As problems 1 and 5 are characterized by a linear 
limit state function and normal probability distribution, the 
application of the methods results in the exact solution of 
the problems. 


From the results obtained in problems 2, 3, 13, 14, and 
15, it can be stated that the FOSM and FORM HLRF 
methods do not converge to stable solutions for highly 
nonlinear limit state functions, with a significant difference 
to the reference (MCBF). In these cases, the FORM 
iHLRF improves the result by decreasing the difference to 
the reliability index reference but still cannot converge to 
the results of the MCBF, MCIS, and SORM methods. The 
high nonlinearity of the functions is the reason for this 
difference between the methods. 


In problem 4, the MCBF method was used with a 
sample of 10º elements and a computational time of 16 
min 35 s and MCIS with 10* simulations with a time of 9 
s, and the FOSM, FORM HLRF, FORM iHLRF, and 
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SORM methods with a time smaller than | s the reliability 
indices of Table 1 were determined. The result of the 
FOSM method is different from the other methods as it 
does not take into account mixed distributions. The MCIS 
showed a better performance for this problem, both in 
computational time and the number of simulations. The 
Monte Carlo, FORM, and SORM simulation methods 
obtained results very close to the reference, with a 
difference of less than 1%. Thus, it can be said that the 
methods MCBF, MCIS, FORM HLRF, FORM iHLRF, 
and SORM adequately represented the values found by the 
authors since the results were very close. 


In problems 6 and 11, only the FOSM method did not 
converge to the result of the other methods because the 
distributions are mixed. Given the weak nonlinearity of the 
limit state function, the FORM HLRF and FORM iHLRF 
methods present identical behavior to the MCBF, MCIS, 
and SORM methods. 


Regarding the results of problem 7, there is a 
difference in the values of the reliability index of the 
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FOSM, FORM HLRF, and FORM iHLRF methods; 
however, the differences found are approximately 1% 
compared to the other methods. This difference presents a 
relatively small divergence. Thus, it can be affirmed that 
this function is weakly nonlinear and that the analyzed 
methods effectively determine the reliability index of this 
function. 


In problems 8 and 10, the MCBF was used with a 
sample with 10º elements and MCIS with 10* elements. 
Despite presenting the same results, the transformation 
methods converge to a result far from the reference 
(MCBF). In these problems, the second-order curvatures 
are zero. A surface of zero curvature is a plane, so FORM 
and SORM got the same results. 


It can be noted that the results obtained in problem 9 
produced very similar results, except for the FOSM 
method. Only the FOSM method did not converge to the 
result of the other transformation methods because it does 
not converge correctly for non-normal distributions and 
still obtained a B greater than the other methods. This lack 
of convergence is because this function is weakly 
nonlinear, and if the other methods had a normal 
probability distribution, they would converge to the same 
result. 


In problem 12, the difference between the reliability 
indices of the Monte Carlo methods compared to the 
FOSM, FORM HLRF, and FORM HLRF is 12%. The 
MCIS once again proved to be more efficient in the 
number of simulations (10* elements) than the MCBF 
(10° elements). 


5.6 
5.2 


4.8 


RELIABILITY INDEX 
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The MCBF method was used with a sample with 107 
elements, MCIS with 10* elements. The analysis of 
problems 16 and 18 was performed considering the MCBF 
method as a reference. SORM and FORM iHLRF 
approximate MCBF and MCIS. There are very marked 
differences between the results of the FOSM and FORM 
HLRF methods with 92% and 78%, respectively, in 
problem 16 to the MCBF reliability index. In problem 18, 
the FOSM and FORM HLRF methods failed to converge 
to the result. The high nonlinearity is the reason for this 
difference. 


The truss presented in problem 17 is a well-known 
benchmark example studied by some researchers. And 
using the MCBF method with a sample with elements, it 
obtained a computational time of 11h 45min 22 s, MCIS 
with 9x/0’ simulations with a time of 15 s, and the FOSM, 
FORM HLRF, FORM iHLRF, and SORM methods with a 
time of less than 5 s. The analysis of the results was 
performed considering the MCBF method as a reference. 
The MCIS method showed once again that it improves the 
performance and efficiency of the MCBF regarding the 
number of simulations and computational time. As 
expected, FOSM did not converge to the benchmark result. 
And the FORM HLRF walked 7% away from the target 
reliability index. 

Figure 15 compares the reliability index of the 18 
problems analyzed. 
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Fig. 15: Comparison of the reliability index of the analyzed examples. 
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It is observed in Figure 15 that the SORM method 
proved to be adequate for the probabilistic analysis and 
returned values equal to or very close to the MCBF and 
MCIS methods in the analyzed problems. SORM is 
calculated by obtaining the reliability index of the FORM 
method. In this way, SORM becomes more advantageous 
because it improves the FORM result by including 
additional information about the curvature of the limit 
state. In the case of problems 8 and 10, the curvatures in 
the failure region are null despite the function being 
nonlinear. 


Also, according to Figure 15, it is observed that the 
MCBF and MCIS methods successfully solved 100% of 
the examples, SORM solved 89%, FOSM, FORM HLRF, 
and FORM iHLRF solved 11%, 44%, and 56%, 
respectively. So, of the transformation methods, SORM 
successfully solved more examples than the FOSM and 
FORM methods. 


Regarding the performance of the FOSM and FORM 
methods, the number of iterations is considered a criterion. 
FOSM and FORM HLRF reached the maximum number 
of iterations (k = 100) in solving problems 11, 13, 14, 15, 
and 16, and FORM iHLRF reached 56 as the maximum 
number of iterations in problem 14. Therefore, there is a 
significant difference based on this performance criterion 
of the FORM iHLRF against the FORM HLRF. 


IV. CONCLUSION 


This work presented a comparative study of structural 
reliability methods (MCBF, MCIS, FOSM, FORM HLRF, 
FORM iHLRF, and SORM), taking into account the 
results of the reliability index. The performance of the 
methods implemented in the Python language was verified 
through eighteen linear, nonlinear and implicit limit state 
function problems. 


The MCBF and MCIS simulation methods worked well 
for all analyzed examples. The MCIS method presented 
the same precision as the MCBF but with greater 
efficiency, converging faster by using a smaller number of 
simulations and with a much lower computational cost in 
analyzing the implicit limit state functions. 


It was observed that the use of the FOSM method is 
quite limited, as it presents results far from the references. 
Thus, it can be said that the technique is not suitable for 
most of the analyzed examples, as it only uses moments of 
first and second order, not taking into account the form of 
the probability distribution. If the basic random variables 
are not normal and the limit state function is linear, the 
imprecision level grows. 
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The FORM HLRF and FORM iHLRF methods did not 
obtain good results for the highly nonlinear functions 
analyzed. Although these methods transform distributions 
that are not normally distributed into equivalent normals, 
they work with a linearization of the failure surface. If the 
function is highly nonlinear around the point of failure, the 
results of this method are not as good as the MCBF and 
MCIS methods. The FORM iHLRF is closer to the 
reference values for the analyzed functions. Therefore, the 
iHLRF algorithm improves the HLRF through a linear 
search for the design point. 


The SORM presented the same or very close results 
concerning the reference (MCBF). The application of this 
method proved to be efficient in sixteen of the eighteen 
problems analyzed and carried out the analysis of second- 
order structural reliability for the functions effectively, 
optimizing in a relevant way the results obtained by the 
FORM iHLRF method. SORM is more beneficial because 
it improves the FORM result by including additional 
information about the limit state curvature at the design 
point. In problems 8 and 10, the second-order curvatures 
are zero. A surface of zero curvature is a plane; FORM and 
SORM gave the same results. 
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